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CN ■ Abstract 
> ■ 

HBT correlations of vr — vr pairs at FAIR/AGS energies are investigated by using the UrQMD 

^ ■ 

00 ' of resonance lifetimes, a soft equation of state and non-Gaussian effects) to understand the HBT 

o ; 

, Rq/Rs puzzle are investigated. Firstly, we find that different treatments of the resonance decay 

^ \ time can not resolve the HBT time-related puzzle, however it can modify the HBT radii at low 

^ ' transverse momenta to some extent to explain the data slightly. Secondly, with a soft equation 

c5 ■ 

■ - - ' of state with momentum dependence, the measured transverse momentum dependent HBT radii 

and Rq/Rs ratio can be described fairly well. Thirdly, non-Gaussian effects are visible in the 
calculated correlation function. Using the Edgeworth expansion, one finds that the non-Gaussian 
effect is strongest in the longitudinal direction and weakest in the sideward direction. 
(Some figures in this article are in colour only in the electronic version) 



transport model and the CRAB analyzing program. Three different possible sources (treatment 
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I. INTRODUCTION 



It is well-known that one can extract information on the space-time dimensions of the 
particle emission source (the region of homogeneity) in heavy-ion collisions (HICs) by using 
the Hanbury-Brown-Twiss interferometry (HBT) techniques. With the ongoing 

advances in detectors and accelerators one is able to explore collision energies for HICs from 
less than ~ 2.5 GeV (SIS/FAIR energy regime), 2.5 - 20 GeV (FAIR/ AGS and SPS) 

up to 20 — 200 GeV (RHIC). Within this broad energy region Quantum Chromodynamics 
(QCD) predicts a transition from a hadron gas to a quark-gluon plasma (QGP), and it is 
expected that this new QGP state of matter exists at least temporarily in the center of HICs. 
During the phase transition (i.e., in the mixed phase), it was proposed that a nontrivial 
transition in the spatio-temporal characteristics of the source exists ^. Unfortunately, so far 
the excitation functions of the HBT parameters have not shown any obvious discontinuities in 
experiments with energies from SIS, AGS, SPS, up to RHIC Instead, several unexpected 
and interesting phenomena occurred, namely, the "E-puzzle" , the "t-puzzle" , and the "non- 



Gaussian" effect 



(for an explanation of these terms, see below). Objectively 
speaking, the existing theoretical investigations are still in-sufficient and further thorough 
investigations are needed. 

In a recent work on the HBT correlation, adopting the Ultra-relativistic Quantum Molec- 



ular Dynamics (UrQMD, v2.2) transport model |9|, [lOl, lUl and the "correlation after-burner" 
(CRAB, v3.0/3) analyzing program 12,[l3|, [l^, the transverse momentum, system-size, cen- 
trality, and rapidity dependence of the HBT parameters of the sources of different identical 
particle pairs (two tts, two Ks, and two As) at AGS, SPS and RHIC energies [?, 



,1151, 



16 



mm 



were investigated. It was found that although the calculations are generally in line with the 
experimental data, discrepancies are not negligible. One of the most puzzling phenomena 
is that the calculated ratio of HBT radii in the outward direction (Rq) and in sideward 
direction (Rs) from central HICs is always larger than that extracted from the data at all 
investigated energies, if the cascade mode is employed, which was named as the HBT time- 
related puzzle ("HBT t-puzzle") 3l- After considering a soft equation of state (EoS) with 
momentum dependence (dubbed as SM-EoS) for formed baryons and a simple Skyrme-like 
(density dependent) potential for the "pre-formed" particles jl^, the HBT radius Rq is 
pushed down and the Rs is pulled up to approach the data so that the "HBT t-puzzle" 
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disappears throughout the whole energy region. Meanwhile, the transverse mass (m^) seal- 



ing^Uch wa. p.edicted in Ret. M, and ha= been p.obed by =eve.a, expenn,ent= recently 
[6|, |20( , can also be much better understood with the help of "pre-formed" particle poten- 
tials in HICs 161]. Therefore, it was concluded that the interaction (here it is embodied with 
potentials) of particles at the early stages of HICs can help to solve the "HBT t-puzzle" . 

However, besides the potential effects, the resonance dynamics and decay may also influ- 
ence the momentum distribution of emitted particles as well as their correlation Since 
this question has been left aside in our previous works, in this paper, we want to complete 
the discussion by investigating the effects of the handling of resonance decay times on the 
HBT quantities by adopting the microscopically transport model UrQMD. The different 
contributions from the resonance decay and from the potential modification are then com- 
pared. In order to have a clearer situation, in this paper we restrict ourselves to the low 
energy region, i.e., the FAIR/AGS energy regime. In this energy region hadronic interac- 
tions dominate the dynamics. String dynamics is negligible and the quarks are still confined 
and the interactions between them do not need to be taken into account during the HICs. A 
further simplification is that the available resonances at these beam energies are mainly the 
A(1232)s, and one can restrict the analysis of resonance modification to this hadron. Due to 
the fact that pions have the largest abundance, they are well suited for the present analysis. 
Negatively charged pions can be easily measured (no contamination from misidentified K"'"s 
and protons) and experimental data are therefore available from most experiments today. 
Thus, we focus the present analysis on negatively charged pions. 

The paper is arranged as follows. In the next section, the UrQMD transport model and 
the CRAB analyzing program are introduced. The different treatments of A decay and 
the effects on pion production are also discussed. In Section 3, firstly, the Gaussian fit- 
ting to the one- and the three-dimensional correlation functions of negatively charged pions 
and the non-Gaussian effect are explored and discussed. Secondly, we show the compar- 
ison of the transverse momentum dependence of HBT radii and Rq/Rs ratios (from the 
Gaussian fitting) between calculations with different treatments of the resonance decay and 
with(without) mean field potentials and experimental data at AGS energies. Finally, in 
section 4, a summary and outlook is given. 



3 



II. MODELS AND TREATMENTS 



A. UrQMD transport model 



The first formal versio n ( verl.O) of the UrQMD transport model was published in the 



end of last century 



21| . Since then, a large number of successful theoretical analyses, 



predictions, and comparisons with data based on this transport model have been accom- 
plished for pp, pA and AA reactions for a large range of beam energies, i.e., from SIS, AGS, 
SPS, up to RHIC [221. 



The UrQMD model is based on analogous principles as the Quantum Molecular Dynamics 
model (QMD) 23, 2^ and the Relativistic Quantum Molecular Dynamics model (RQMD) 



251]. Similar to QMD, hadrons are represented by Gaussian wave packets in phase space. 



and the phase space of hadron i is propagated according to Hamilton's equation of motion: 



OH 



dH 



Here r and p are the coordinate and momentum of hadron i. The Hamiltonian H consists 
of the kinetic energy T and the effective two-body interaction potential energy V, 



H = T + V, 



(2) 



and 



T = Y.iE. - m,) = 5:(v'mf + - m.). (3) 

In the standard version of UrQMD model [o], Q, 21], the potential energies include the 
two-body and three-body (which can be approximately written in the form of two-body 
interaction) Skyrme- (also called as the density dependent terms), Yukawa-, Coulomb-, and 
Pauli-terms as a base. 



V = , + VjUl + VYuk + Vcou + Vpau. 



(4) 



sky ' sky 

And the single particle potential follows from U = 6V/6f, where / is the phase space 
distribution function which reads as 



/(r,p) = E/^(r,P) = E 



(5) 
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here L is the width parameter of the wave packet. Recently, in order to be more successfully 
applied into the intermediate energy region (0.1 ^ Eb ^ 2A GeV), the UrQMD has contained 
more potential terms 26|], those are, 1) the density-dependent symmetry potential, which 



is essential for isospin-asymmetric reactions at intermediate and low energies, and 2) the 
momentum-dependent term {2?! 



Umd = tmdln^[l + amdiPi " Pj)^]Pi/PO 



(6) 



where tmd and amd are parameters, po is the normal density, pi is the density of the baryon 



pi 



p(ri)pdr= / p(ri)^p(rj)dr 



(47rL)3/' 



(ri-n)^ 



(7) 



J ^ ' J 

With these updates, some sensitive probes of the (density dependent) symmetry potential 



have been proposed 



26 



28 



29 



30[ |. Further, it was also found that the experimental pion 



and proton directed and elliptic flows from HICs with beam energies from ~ lOOA MeV to 
2A GeV can be well described 

At higher beam energies, the Yukawa-, Pauli-, and symmetry- potentials of baryons 
becomes negligible, while the Skyrme- and the momentum-dependent part of potentials still 
influence the whole dynamical process of HICs. With the help of a covariant prescription of 
mean field from RQMD/S {32!, the effects of the mean field with momentum dependence on 
collective flows from HICs at 2-1 58A GeV energies were studied by a Jet AA Microscopic 
Transportation Model (JAM) and it has been found that the momentum dependence in the 
nuclear mean field is important for the understanding of the proton collective flows at AGS 
and even at SPS energies 33|. In this work, we choose the same SM-EoS as that in Ref. jsS], 
where the momentum dependent term reads as 



md 



E 



''md 



/(ri,Pj 



[(Pi-Pi)/aL] 



(8) 



k=l,2 Po 

where the Umd consists of two parts with separate parameters in order to fit better the real 
part of the optical potential 3J]. We found that the momentum dependent term is also 
important to explain the transverse momentum dependence of the HBT parameters at AGS 
energies 
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35l |. Furthermore, as in Ref. j33j, the relativistic effects on the relative distance 



Tij = Ti — Tj and the relative momentum pij = Pi — pj employed in the two-body potentials 
(Lorentz transformation) are considered as follows: 



for A I 


1=1 /r 
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FIG. 1: Lifetime r of the A(1232) resonance defined as the inverse of resonance width. 1/2, 1, and 
2 times of A(1232) width at the pole mass are selected as well as a mass dependent one (r(M)). 



^ = 4+%(^v■^3^J)^ (9) 



p^, = Pl-iE.~E,f+^U%^r- (10) 

In Eqs. M and [10] the velocity factor jSij and the corresponding 7-factor of i and j particles 



are defined as Pij = (pj + pj)/{Ei + Ej) and = 1/y^l — [3fy 

A covariance-related reduction factor for potentials in the Hamiltonian, rrij/Ej, was in- 
troduced in the simplified version of RQMD model 32] and adopted in this work as well. 

The collision term of the UrQMD model treats 55 different baryon species (including 
nucleon, A, A, S, S, and Q and their resonances with masses up to 2.25 GeV) and 32 
different meson species (including light unflavored and strange mesons and their resonances 
with masses up to 2.0 GeV) as tabulated in the PDG [36]. Through baryon-antibaryon 
symmetry the respective antibaryons are included. The isospin is explicitly treated as well 
jy, [l^- For hadronic continuum excitations a string model is used. Starting from the version 
2.0, the PYTHIA is incorporated into UrQMD in order to investigate the jet production and 
fragmentation at RHIC energies 11|. 
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B. Treatments of the A resonance decay in UrQMD 



The frequently used approach to resonance hfetimes in most transport calculations is the 
application of r = l/F in conjunction with a Monte-Carlo sampling of the exponential decay 
law. The decay width F of a resonance is usually either taken to be constant (Fq) or mass 
dependent (F(M)) 9,[l^. Fig. [T]shows the lifetimes of the A(1232) resonance as a function 
of its invariant mass. Considering that the spectral functions might be modified (width 
and the position of the mass pole) by the hot and dense nuclear medium - e.g., recently a 
precision measurement of low-mass muon pairs in 158 A GeV Indium-Indium collisions at the 
CERN-SPS was reported, and the associated space-time averaged p-meson spectral function 
shows a strong broadening, but essentially no shift in mass [54] - it is desirable to explore 
the effect of the broadening/ narrowing of resonances on pion freeze-out. Therefore, two 
other mass-independent lifetimes for A(1232) are adopted which correspond to the widths 
l/2Fo and 2Fo, respectively (shown also in Fig. [1]). The lifetimes of other resonances are not 
altered in the calculations. 

It should be noted that the treatment of the resonance lifetimes and widths in a medium is 
a long standing problem and much argued in the theory community. In fact, the widely used 
prescription for resonance lifetimes (especially for broad resonances) in transport models 
are still under debate. Especially, the recent definition of the lifetime by the derivative 
of phase shift with respect to the center-of-mass energy is thought to be more reasonable 



38, 



39 



40| . However, there still exist difficulties to adopt the "new" lifetime consistently 



in transport models. In addition, as a result of the scattering, the forward going part of 
the wave packet can suffer a different time delay from the scattered one, which is not trivial 
to implement, and, it is still not understood how to tackle the issue of negative time delay 
in the real transport model calculations j^, Q. Due to this ongoing debate, we stick to 
the above mentioned schematic change by factors of 2 to see if any effect emerges from this 
modification at all. 

Fig. [2] (upper plots) depicts transverse mass spectra of charged pions at midrapidity 
{\Ycm\ < 0.05, Ycm = |log( ^'""^p|| ), Ecm and p\\ are the energy and longitudinal momentum 
of the pion in the center-of-mass system) for central (< 5% cr-r) Au+Au collisions at 2A 
GeV. The lower plot shows the ratio between tt^ and 7r+ transverse mass spectra. The data 
are taken from Ref. 4l|- The lines correspond to calculations with and without potentials 
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Cas, r(M) 
SM-EoS, r(M) 



0.8 

m,-m^ (GeV/c ) 



FIG. 2: Upper plots: Charged pion transverse mass spectra at midrapidity (|ycm| < 0.05) for 
central (< 5% ctt) Au+Au collisions at 2A GeV. Lower left plot: The ratio between vr" and 7r+ 
transverse mass spectra. The calculations with various lifetimes in cascade mode ("Cas") as well 



4l|. The horizontal lines 



as in potential mode ("SM-EoS") are compared with the E895 data 
represent the unity and the value (5A^^ + N P) / {5P'^ + NP), separately (see text). 

(the "SM-EoS" and the "Cas" modes), as well as calculations with mass-dependent lifetimes 
for all resonances (with respect to "r(M)") and with three mass-independent lifetimes for 
A(1232) resonance (with respect to 1/2T^, T^, and 2r^, respectively) while other particles 
are not changed. In line with the previous UrQMD calculations [11!], the vr transverse mass 
spectrum can be well described by the cascade calculation, except for a slight deviation 
from the data at large transverse masses. It is found that a broader constant width of the 
A(1232) resonance makes the spectra a little steeper since the A resonances decay earlier. 
The mass dependence of the resonance lifetime only influences slightly the spectra in the 
low transverse mass region, which is mainly due that the A resonances with their relatively 
small masses have long lifetimes, as a consequence, have a higher probability to be absorbed 
through the A^A —>■ NN detailed balance process. However, it is interesting to see in the 
lower plot of Fig. [2] that the ratio between vr" and vr"*" spectra is reproduced fairly well with 
the "SM-EoS" mode, while there is no transverse mass dependence in the cascade modes. 
The value {5N^ + NP)/{5P^ + NP) ("A^" and "P" represent the initial neutron and proton 
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numbers of HICs) of the horizontal hne is obtained if only the single pion production via 
A resonances in nucleon-nucleon collisions is considered. However, the effects of such as 
potentials, reabsorptions, and rescatterings of pions. As, and other resonances will change 
the value visibly. The strong transverse mass dependence of the 7r~/7r+ ratio observed in 



this plot is mainly due to the Cou 
investigated before 



26 



27 



30 



42 



omb interaction between charged particles which has been 



43|. 



C. CRAB analyzing program and the fitting process 

To calculate the two-particle correlator, the CRAB program is adopted which is 
based on the formula: 

^ ^ J d^Xid'^Xjg{xi,pi)g{xj, pj)\(f){r' , q^)\'^ 
J d^Xig{xi,Pi)J d^Xjg{xj,Pj) 

Here g{xi,pi) is an effective probability for emitting a particle i with 4-momentum pi = 
(i?j,Pj) from the space-time point Xi = (tj,rj). (/)(r', q') is the relative wave function with 
r' being the relative position in the pair's rest frame, q = Pi — Pj and k = (pj + Pj)/2 are 
the relative momentum and the average momentum of the two particles i and j. Due to 
the quantum statistics, the correlator is larger than unit at small q for bosons, and in the 
absence of strong and Coulomb final state interactions, the wave function of an identical 
pair of bosons reads as 

|</)(r',q')P = l + cos(2q'-r'). (12) 

And the two-particle correlator and the "region of homogeneity" can be directly related by 
a Fourier transformation. However, in HICs, the Coulomb and strong interactions might 
distort the wave function and the direct relationship between the correlator and the region 
of homogeneity disappears. As a standard method, one can fit the correlator as a three- 
dimensional (3D) Gaussian form under various reference frames, in which the longitudinal 
comoving system (LCMS) (or called as the "Out-Side-Long" system) is more frequently 
adopted in recent years. We use it in the present work as well. The corresponding 3D 
Gaussian correlation function can be expressed as 

C{qo, qs, Ql) = K[l + Aexp(-i?ig2 _ ^2^^2 _ ^2^2 _ 2i?2 ^go^L)]- (13) 
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In Eq. [13] the K is the overall normalization factor, the qx and Rx are the components 
of the pair relative momentum and homogeneity length (HBT radius) in the x direction, 
respectively. The A parameter is commonly called the incoherence factor and lies between 
(complete coherence) and 1 (complete incoherence) for bosons in realistic HICs. Because 
the parameter A might be influenced by many additional factors, such as contamination, 
long-lived resonances, or the details of the Coulomb modification in the FSI, we regard it 
as a free parameter. The represents the cross-term and plays a role at large rapidity 

0,0, Ci- 



However, due to the resonance decay and the space-time correlation [45|, l46|, |47|] , the emis- 
sion function has been found to deviate from a Gaussian form (the non-Gaussian effect). 
In order to gain clearer information about the source, both fitting the correlator to a (3- 
dimensional) Gaussian form with higher orders of a harmonics and the (3-dimensional) imag- 

le high orders of harmonic fitting by the Edgeworth 



ing technique [48|, |49| are available. T 



expansion was proposed by Csorgo [50| and used in experiments [511, which is expressed as 



C^(ginv) = i^inv[l + Ainv exp(-i?LgL)] [1 + E ^77=^^n(^mvginv)]. (M) 

Here the subscription "inv" means that the one-dimensional (ID) correlation function is 



constructed in the invariant quantity ginv = y {(fY — (q)^ while -Rmv is the corresponding 
ID radius, n is the order parameter, and the if„ reads as 

Jn 

H4z) = i-ire^' —e'^\ (15) 

To calculate the HBT two-pion correlation, firstly, we select central collisions at AGS 
energies in the UrQMD model: Au-^Au at Eb = 2, 4, 6, and 8A GeV (< 11% of the total 
cross section ctt), with a rapidity cut \Ycm\ < 0.5. For each case about 40 thousand events are 
calculated. All particles with their phase space coordinates at their freeze-out times are put 
into the CRAB analyzing program. Only the negatively charged pions are considered during 
the analyzing process in this work (for each analysis, one billion pion pairs are considered . 



We found that the Coulomb effect in FSI on the HBT radii of the pion source is small 
so that we do not consider it in the calculations. 
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FIG. 3: Projections of the correlation function for negatively charged pions in longitudinal (top 
plots), outward (middle) and sideward (bottom) directions from central Au+Au collisions at -E^ = 2 
A GeV. The results with "Cas" (left plots) and "SM-EoS" modes (right) are shown with points. 
The mass-dependent lifetime of resonances is employed in the calculations. ID Gaussian fit to each 
projection is shown with lines as well. 

III. HBT RESULTS 

A. Gaussian fitting 

In Fig. [3]we show the projections of the calculated correlation function of the pion source 
in longitudinal (top plots), outward (middle) and sideward (bottom) directions without 
( "Cas" , squares in left plots) and with ( "SM-EoS" , circles in right plots) potentials. The mass 
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FIG. 4: ID correlation function for vr^vr^ pairs. Fittings with the Edgeworth expansion up to the 
eighth order on it are also shown. The result with cascade mode (left plot) is compared to that 
with the SM-EoS mode (right). The fitting results of the radii Rirw are listed in the plots as well. 

dependent lifetime of resonances is adopted. Correspondingly, ID Gaussian fits {C{qx) = 
Kx[l + Xx Qw{~RIQx)] where "x" stands for "L", "O", or "S"') to them are also shown with 
dashed and dotted lines, separately (the values of the HBT radius in each projection and the 
corresponding reduced least square fit S/dof are also given in each plot). The fitting result 
shows clearly that the non-Gaussian effect of the correlation function exists especially in 
the longitudinal direction, while in other two directions it is largely reduced. The Gaussian 
fitting in the sideward direction is even better than that in the outward direction. We notice 
that in the data analysis by both the HBT fitting and the imaging method the similar 
non-Gaussian effect has been observed as well at all AGS, SPS, and RHIC beam energies 



51 



52 



531]. The non-Gaussian character originates mainly from both the decay of long- and 



intermediate-lived resonances, such as r]', uj,_ K* mesons, and A, A*, E* baryons, and the 
space-time correlation (or fiow) effect 45 



46l . |47| . At AGS energies, the contributions from 
meson and hyperon resonances are limited, but the A resonance decay plays an important 
role. The origin of the non-Gaussian effect shown in Fig. [3] might be clear by comparing the 
HBT results with and without the mean field potential. If comparing the left and the right 
plots of Fig. [3] roughly, one may find that the Gaussian fitting becomes even better with the 
consideration of the mean field potential. 
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To quantitatively analyze the non-Gaussian effect, the Edgeworth expansions up to the 
8th order on the ID correlation function are shown in Fig. |H It is seen that the -Rinv 
increases with the order of the expansion on the correlation function. In the cascade mode, 
the difference of the i?inv values between without expansion and with expansion up to the 
gth Qj-f^Qj- reaches 1.18 fm, but in the SM-EoS mode the difference is reduced to 0.83 fm. 
Similarly, if one takes the same expansions into the long-out-side projections, it is found 
that, 1) in the cascade mode, the radii in longitudinal, outward, and sideward directions 
change 2.2, 0.6, and 0.52 fm, respectively. 2) While in the SM-EoS mode, the changes 
are 1.45, 0.51, and 0.38 fm, respectively. Therefore, the non-Gaussian effect is reduced in 
all directions of the correlation after considering the mean field potential. Therefore, it 
is the long- and intermediate-lived resonance decay that mainly contributes to the visible 
non-Gaussian phenomenon but not the mean field. 

B. Transverse momentum dependence of the HBT radii 

Fig. Ogives the transverse momentum kx (kr = (Pit + P2t)/2) dependence of the 
HBT-radii Rl, Ro, Rs, and the Rq/Rs ratio in central Au+Au collisions at Eb = 2A GeV. 
Experimental data are taken from Ref . {35! • The cascade calculations with mass-independent 
lifetimes of resonances are shown with different widths of the A(1232) resonance (the widths 
of other resonances are not changed). With the default width (Fq), we find that the radii Rl 
and Rs are somewhat smaller than the data at small kx. With a narrower A(1232) width 
(hence a longer lifetime) the HBT radii at large kx are increased slightly, and vice versa. 
This modification can not explain the experimental Rq and Rs data, which can also be 
clearly seen from their Rq/Rs ratios in the right-bottom plot: the variation of the A(1232) 
width does not help to obtain a small value of the Rq / Rs ratio as observed in the data. 

Comparing the left and the right plots of Fig. [3] more carefully one can find that the mean 
field potential enhances the incoherence in the reaction plane (x-z plane, correspondingly, 
the outward-longitudinal plane) but not in the sideward direction. Fig. [6] illustrates the kr 
dependent HBT radii Rl, Rq, Rs, and Rq/Rs ratio (from top to bottom plots) of the pion 
source for the energies Ef, = 2, 4, 6, and 8A GeV (from left to right plots). With the mass 
dependent lifetimes for all resonances in the cascade mode ("Cas,r(M)"), the resonances 
with their small invariant masses decay later and hence expand the source when compared 
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FIG. 5: Transverse momentum fcy dependence of the HBT radii Rl, Rq, Rs^ and the Rq/Rs ratio 
in central Au+Au cohisions at Ei, = 2A GeV. The cascade mode is chosen. The mass-independent 
Ufetime of resonances is considered, while the width of A(1232) resonance is varied from l/SFn 



^1 



To, to 2ro (the widths of other resonances are not changed). Experimental data are shown 
with scattered stars. 

with the mass-independent lifetime of resonances ("Cas,ro"). Considering the fact that 
the resonances with small invariant masses tend to produce pions with small momenta, an 
increase of the HBT radii at small and Ef, is understandable. At large kj- as well as at 
high beam energies, this effect is reduced since less resonances with small invariant masses 
contribute. It is interesting to see that the result with a mass dependent treatment of 
resonance lifetimes can match the HBT data better in the AGS energy region. However, it 
is also seen that the mass dependence will increase the HBT radii with almost equal power in 
outward and sideward directions and the ratio between Rq and Rs values remains basically 
unchanged. With the help of the SM-EoS, it is interesting to see that the Rq at large kx 
are driven down while the Rs at small /ct are pulled up so that the /ct dependence of the 
Ro/Rs ratio can be much better described, in line with Ref. [15|- This finding supports 
that it is the mean field which leads to a stronger phase-space correlation that results in a 
much better explanation of the HBT time-related puzzle (l6l |. 
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FIG. 6: dependence of HBT radii Rl, Rq, and Rs as well as the Rq/Rs ratio of negatively 
charged pion source (from top to bottom plots) for central Au+Au collisions at = 2, 4, 6, and 
8A GeV (from left to right plots). In calculations, 1) the"Cas,ro" represents the cascade mode 
with a mass-independent lifetime of resonances. 2) "Cas,r(M)" means the cascade mode with 
a mass-dependent lifetime of resonances. 3) While the "SM-EoS,r(M)" represents the potential 
mode with a mass-dependent lifetime of resonances. The experimental data are taken from jsS]. 



IV. SUMMARY AND OUTLOOK 



To summarize, the HBT correlation of vr^ — pairs at AGS energies were investigated 
by using the UrQMD transport model and the CRAB analyzing program. We found that 
a narrower (wider) width of A(1232) resonance results in larger (smaller) HBT radii in all 
directions at large k^, which can not explain the small experimental ratio between Rq and 
Rs- Although a mass dependent lifetime of the resonances can not resolve this problem as 
well, but it pulls up the HBT ratii at small kr and hence slightly improves the HBT radii 
of pions. 

We observed that the kr dependent HBT radii and Rq / Rs ratio in the AGS energy region 
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can be described fairly well with a soft equation of state with momentum dependence. This 
supports the idea that the interaction of particles in the early stage of the reaction (leading 
to stronger correlation) is the key to solve the HBT time- related puzzle. 

Non-Gaussian effects are visible (although weak) in the correlation function and might 
bring large uncertainties (on the order of ±1 fm) if the correlation function is fitted only by 
a Gaussian form. To investigate the non-Gaussian effect, the Edgeworth expansion was used 
during the fitting process. It was found that the non-Gaussian effect is strongest in the longi- 
tudinal direction and weakest in the sideward direction. The decay of the intermediate- and 
long-lived resonances was found as the main contribution to the non-Gaussian phenomenon, 
while the mean field potential did not increase this effect. 

At higher beam energies, such as SPS and RHIC, the non-Gaussian effects as well as 
the broadening of vector mesons (in connection to the NA60 results js^) and its influence 
on TT+TT^ correlation deserve further investigations. In addition, other, besides pions, (non- 
) identical particle correlations deserve increased attention as they might provide snapshots of 
the source geometry at different times of the reaction. These studies are currently underway 
and will be addressed in a forthcoming paper. 
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